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Abstract. We present a geometric characterization of the non- 
linear smooth functions V : R —> K for which the origin is a global 
isochronous center for the scalar equation x = —V'(x). We revisit 
Stillinger and Dorignac isochronous potentials V and show a new 
simple explicit family. Implicit examples are easily produced. 
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^ ' Dedicated to Jorge Lewowicz for his 75 th birthday 



1. Introduction 



O ' Consider the scalar differential equation x = —g(x) under the mini- 

mal assumptions that g is continuous, g'(0) exists and it is > 0, and that 
xg(x) > whenever x^O. Let V be a primitive of g, with V(0) = 0. 
As it is well-known, the Cauchy problems have unique solution, thanks 
to the energy first integral x 2 /2 + V(x), all solutions with small enough 
^ . energy are periodic, and their orbits in the x, x plane encircle the origin 

(0, 0), which is called a center for the system. If moreover V(x) — > +oo 
as x — > ±oo, then all solutions are periodic, and the origin (0, 0) is a 
global center. When all orbits around a (global) center have the same 
t*** . period the (global) center is called isochronous. 

^ j After the 1937 paper [8] by Koukles and Piskounov, a construc- 

tive characterizations of local isochronous centers was proposed by 
Urabe [T3| [H] in 1961-62, in terms of a differential equation that must 
be satisfied by V, and that contains an arbitrary odd function. Urabe's 
method has been cited and developed by many authors, for example 
CN ] recently by Chouikha, see Strelcyn [T2] and the reference therein. 

In 1989 the paper [TB] gave a completely different characterization of 
local isochronous potentials (Theorem 12.31 below) . It is more construc- 
^ | tive than Urabe's method, because we can start with any involution h 

(Definition 12. II below) . plug it into the simple formula c(x — h(x)) 2 and 
the result is an isochronous potential. The involution condition has a 
very clear geometric meaning: the graph of h must be symmetric with 
respect to the diagonal. The involution criterion was independently 
rediscovered by Cima, Manosas and Villadelprat [5] in 1999 among 
several other results on more general Hamiltonian systems. 

In the last few years the literature has seen some resurgent interest in 
general isochronous systems (Calogero's book [3] is a prime example), 
and particularly in characterizing global isochronous potentials and in 
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exhibiting elementary explicit examples, beside the trivial linear oscil- 
lator. To our knowledge, the earliest global analytic example was given 
by Stillinger and Stillinger [TT] in 1989 in a journal addressed to the 
chemical community. A subclass of Stillingers' example was rediscov- 
ered by Bolotin and MacKay [2] in 2003, together with necessary and 
sufficient conditions for a global center. Dorignac [7] in 2005 was the 
first to cite the Stillingers, produced new examples and noticed that we 
can introduce a rescaling parameter. Among the most recent papers 
we may cite the ones by Asorey et al. [1] and by Mamode [9]. 

Two papers that concern us most closely appeared in 2011. The 
first [T7| gave a hopefully more accessible account of the involution 
method, and then it connected it to some Hamiltonian systems in di- 
mension 4 which are interesting for Lyapunov stability, also finding lo- 
cal isochrony examples by means of that dynamics (see also Section H] 
below). The second paper [6], by Cima, Gasull and Manosas, included 
the construction of an explicit global family of nonlinear centers, which 
was shown to be isochronous by the technique of involutions too. It 
was this paper which first drew our attention to the problem of global 
examples. 

The purpose of this note is to further demonstrate how the involution 
method is very well suited for the task of finding plenty of explicit global 
isochronous centers. The original paper [16] suggested that involutions 
can be found by rotating the graph of an even function, but it kept the 
discussion local and did not pursue examples. Here we describe another 
geometric idea which seems preferable for global purposes: we look for 
involutions as (subsets of) the zero set of two- variable functions f(x, y) 
which are symmetric, i.e., f(x,y) = f(y,x) (Theorem 13.11) . Two new 
examples will be readily obtained this way in Section El 

The idea of using symmetric f(x,y) to obtain involutions is elemen- 
tary, and we are thankful to Armengol Gasull for pointing out to us 
that something similar, although with applications unrelated to ours, 
was described in the 2002 paper [15], which in turn draws from the 
1967 [ID]. 

In Section |4] we show that Stillingers' examples can be obtained by 
choosing / to be a symmetric quadratic polynomial, whose zero set is 
a hyperbola. It also shows that an example from [T7] contains Still- 
ingers'. Section [5] finally finds a symmetric function / that yields one 
of Dorignac's examples that we revisit because of its remarkable sim- 
plicity. 

2. Constructing isochronous centers 

The concept of involution is crucial in our approach. Basically, a 
function h : A — > A is an involution if hoh is the identity mapping on A. 
In this paper we will use the term "involution" to mean something more 
specific, as follows: 
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Definition 2.1. A C 1 diffeomorphism h of an open interval J C R 
onto itse// is called an involution if 

(2.1) /r 1 = h, G J, h(0) = 0, h'{0) = -1. 

Geometrically speaking, the condition = /i means that the graph 
of h is symmetric with respect to the diagonal; indeed (x,h(x)) has 
(h(x) , X^j clS symmetric point and this coincides with the point (h(x), 
h(h(x))) of the graph. The condition h(0) = means that the graph 
intersects the diagonal at the origin. For example, start with the hy- 
perbola yx = l, which is symmetric with respect to the diagonal. If we 
translate its point (1, 1) to the origin we get (y + l)(x + 1) = 1, which 
can be solved for y as y = —x/(l + x). If we finally take the branch 
that goes through the origin we arrive at the following involution 

(2.2) h(x) = — , x G J = (-1, +oo) . 

1 + x 

A homothety of the plane transforms an involution into another in- 
volution: 

Remark 2.2. Let a G R \ {0} and h be an involution on (b,c), then 
h(x) = h(ax)/a is an involution on (b/a,c/a) if a > 0, on (c/a,b/a) 
otherwise. 

In particular h x = h, indeed h(h(x)) = h(ah(ax)/a)/a = x. In 
this way f)2.2p gives the following 1-parameter family of involutions 

{(-1/a, +oo), a>0 
(— oo, +oo), a = 
(— oo, —1/a), a < 

These are the only involutions that are rational functions of x, as proved 
in 0. 

Theorem 2.3. Let h : J — )■ J be an involution, u > 0, and define 

(2.4) V(x) = y (x-h(x))\ xeJ. 

Then the origin is an isochronous center for x = —g(x), where g(x) = 
V'(x), namely all orbits which intersect the J interval of the x-axis 
in the x,x-plane, are periodic and have the same period 2tt/u. Vice 
versa, let g be continuous on a neighbourhood of G R, g(0) = 0, 
suppose there exists g'(0) > 0, and the origin is an isochronous center 
for x = —g(x), then there exist an open interval J, G J, which is a 
subset of the domain of g, and an involution h : J — >■ J such that ( 12. 4p 
holds with V(x) = J* g(s)ds and oo = a/o'(0). 

The proof is included in the proof of Proposition 1 in [16], by one 
of the present authors, as a particular case. Formula ( 12. 4p corresponds 
to formula (6.2) in the paper [16]. A detailed proof can be also found 
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in the recent [T7]; see Theorem 2.1 and Corollary 2.2 in [T7]. This 
last paper also contains the following local necessary conditions for 
isochrony: 

(2.5) V (U) - 3y//(Q) V ([)) — ^ gy//(0)3 , 

which can be deduced by taking successive derivatives of the involution 
relation h(h(x)) = x at x = 0. 

Inserting the involution (12. 3p into formula (12 .4p we obtain the fol- 
lowing isochronous potential 

2 «, 2 f(-l/a,+oo), a>0 

(2.6) V{x) = ^x 2 (^^) , xe J= j(-oo,+oo), a = 

I (— oo, —1/a), a < 0. 

This is the only rational potential for which there is isochrony, as proved 
in [1] and in [TJ. Of course, it is not defined globally on R, except in 
the trivial case a = 0. 

Formula (12 .4p implies the following global inequality: 

u 2 

(2.7) Vfa) > —x 2 , xeJ. 

8 

Simply observe that x and h(x) have opposite signs, so that x—h(x) > x 
when x > 0, and x — /i(x) < x when x < 0. The inequality ( 12. 7p is 
meaningful not for small x (because V(x) ~ u 2 x 2 /2 as x — > 0) but 
rather as x — >■ ±oo. In particular, global isochronous potentials are 
at least quadratic at infinity. If V(x) grows more than quadratically 
as x — ¥ +oo, it must grow not more than quadratically as x — ¥ — oo, 
and the other way round. 

3. Global centers 

Our approach to centers on the whole M. 2 for x, = —g(x) is based on 
the following geometric characterization of global involutions in terms 
of zero sets of symmetric functions of two variables: 

Theorem 3.1. Let f : Q —¥ R be a C 1 function on the open set flCR 2 
such that: (0, 0) G Q, f(0, 0) = 0, and 

(3.1) fay) e SI => (y,x)eSl, f(y,x)=ffay). 

Let T be the connected component o// _1 (0) that contains the origin. 
Suppose that d 2 f(x,y) ^ for all (x,y) G T, and that T project itself 
onto the whole x-axis, namely, for each x G K we have (x, y) G V for 
some y G R. Then T is the graph of a global involution h that is defined 
on the whole R. All global involutions can be obtained this way. 

Proof. By the implicit function theorem, T is the graph of a C l function 
h with h(0) = 0. Moreover, h is defined on the whole R since T projects 
itself onto the whole From ( 13. ip we have d 1 f(x,y) = d 2 f(y,x) 
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Figure 1. Global involution denned by (13. 2p . 



for (x, y) G Q so d\f never vanishes on Y and has the same sign as d 2 f ■ 
We deduce that 

,, x dif(x,h(x)) 

h & = - » 7f \) i\ < °' x e Rj 

d 2 f(x,h(x)) 

and in particular /i'(0) = —1. Finally, f(h(y),y) = f(y,h(y)) = 
shows that = h and ft, is a global involution. 

Conversely, let h : R — > R be a global involution and define /(x, y) : = 
x + y — h(x) — h(y). This is a C 1 function on R 2 , /(0, 0) = 0, and 
/(y, x) = /(x, y). We have d 2 f(x, y) = 1 — /i'(y) > for all (x, y) € M 2 . 
The graph of h is a subset of the set T in the statement of the theorem, 
so T projects itself onto the x-axis (and it actually coincides with the 
graph of h). □ 

The theorem can be applied to the / of the form f(x,y) = y?(x) + 
<f(y), where (p is a diffeomorphism of R onto itself, and <p(0) = 0. 
For a non trivial example let us consider the following function of two 
variables 

(3.2) /( Xj3/ ) : = a; _ a ;2 +a .5 +3/ _ 3/ 2 +3/ 5 j ( x ,y) e m 2 . 

All conditions in Theorem 13.11 are satisfied since d%f(x,y) = 1 — 2y + 
5y 4 for all y G R and the function x i— > x — x 2 + x 5 is a diffeo- 
morphism l->lso the zero set of / projects itself onto the whole R. 
So a global (implicit) involution h is defined. It is a non elementary 
algebraic function. In this example the zero set is connected and it 
coincides with T. The graph of h is depicted in Figured! 

To get an example where there are more connected components, and 
d 2 f(x,y) = at some (x, y) T, we can simply modify / as 

f(x,y) = (x-x 2 + x 5 + y-y 2 + y 5 )((x-l) 2 + (y-l) 2 -l), 
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Figure 2. Global involution ( 13. 4p . 

which yields the same global involution (notice that c^/ (0, 1) = 0). 

Sometimes the equation f(x,y) = can even be solved explicitly. 
For instance for the family 

(3.3) f(x, y):=x + y + p(e x + e y -2)=0, (x, y) G R 2 , p > 0. 

Again, all conditions in Theorem 13.11 are satisfied since &2f(x,y) = 
1 + pe y for all y G R and the function x i— > x + pe x is a diffeomorphism 
R — > R so the zero set of / projects itself on the whole R. In this 
case the global involution can be written in terms of the Lambert real 
function Wo : [— 1/e, +oo) — > [— l,+oo) (the principal branch of the 
product logarithm so Wo(xe x ) = x), as 

x^-x + p{2- e x ) - W (p e ~ x+p{2 - eX) ) . 
By Remark 12.21 we finally obtain the two real parameters global family 

(3.4) h(x) = - x + £(2- e ax ) --W (p e —+p(^ ax )\ 

a a ' 

x e R,p > 0,a ^ 0. 

In Figure [2] we can see the graph of the involution h for a > 0. The 
growth of h(x) is exponential as x — » +oo and logarithmic as x 7 
— oo. Hence the corresponding isochronous potentials V(x) will be 
exponential asx-)> +oo and quadratic as x — > — oo. 

4. Stillinger centers from hyperbolas 

Let us consider the following symmetric quadratic function with two 
real parameters 

(4.1) f(x, y) = \{x 2 + y 2 ) + 2(A + 2a 2 )xy + 4a(A + a 2 )(x + y), 

(x,y)eR 2 , A>0, a^0. 
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The equation f(x,y) = gives a hyperbola with center (—a, —a), ver- 
tices (0, 0), (—2a, —2a), and asymptotes 



(4.2) X(y + a) = (-(A + 2a 2 ) ± 2aV\ + a 2 ) (x + a). 



Both values — (A + 2a 2 ) ± 2a VA + a 2 < 0, so each branch of the hy- 
perbola projects itself on the whole x-axis. If a > the upper branch 
intersects the origin, the lower one if a < (in both cases for x — > — oo 
the asymptote of the branch through the origin is the one with the + 
sign in (14. 2p and for x — > +oo the one with the — sign). 

All conditions in Theorem 13. II hold, in particular f(0, 0) = 0, f(x, y) 
= f(y,x) for all (x, y) G M 2 , and we can check that d2f(x,y) ^ for 
all points (x,y) such that f(x,y) = 0. By Theorem 13.11 the equation 
f(x,y) = defines a global involution h : R — > M. In the present case 
we can even write down explicitly 



(4.3) h(x) := -a- ( (A + 2a 2 )(o 4 .r,H 



2aVA + a 2 a/A + (a + a;) 2 ) . 



x G E, A > 0, 



where a G M, since a = gives the trivial h(x) = —x. 

Replacing the involution into formula (12. 4p we obtain the following 
isochronous potential 

(4.4) V(x) = u 2 ^ A a 2 + (a + x) ((A + 2a 2 )(a + x) + 

- 2aVA + a 2 a/A + (a + x) 2 ^ , xGM,w^0,A>0. 

So the origin is a global center for i = —g(x) with g(x) = V'(x), i.e., 

. 2 (A + a 2 )l/. . A + 2a 2 A + 2(a + x) 2 

(4.5) g(x)=Lo- — — — (a + x) — ■ — a 
1 ^ 1 ^ A 2 1/ ; a/AT^ A/A-Ka + x)^ 

The function g is linear if and only if a = 0, in which case g(x) = u 2 x. 

Remark 4.1. If we introduce the new parameters 

(4.6) A = ^T"' a = 7^' ee(-i,i),/3>o, 

t/ien (14. 4p ta&es Stilling er's form (see [11] formula (4)) 



(4.7) 



a; 2 



2(1 -e 



(x + J|=(\-^1 + 2£a/^x + /3x 2 ^ . 



A totally different technique was used in [T7] to find some explicit 
examples of local isochronous centers. We are going to show that one 
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of them can be rewritten as formula (14.51) after a suitable reparameter- 
ization. Indeed, consider the g defined in formula (5.12) in [T7], that 
is, 

/ x / \ 2cw 2 f /l2 ,„ . , 6 2 - 4c - 2(6 + 2x) 
(4.8) g(x) = — — (6 2 + 4c) (6 + 2x) + b ■ 



(b 2 -Ac) 2 y ^/l + x(b + x)/c 

Here u > 0, b, c G R, c ^ 0, b 2 — 4c 7^ 0. Each of these functions g 
in (14 .8p gives a global center on the whole R 2 whenever it is defined 
on the whole R, that is, if and only if 1 + x{b + x)/c > 0, so that its 
square root in the denominator in (14. 8ft exists for all x G R and never 
vanishes. This is equivalent to b 2 — 4c < 0. By introducing the new 
parameters A := 4c — b 2 , a := 6/2 we have that the condition A > 
is necessary and sufficient for g to be defined on the whole R. Notice 
that this implies c > 0, in particular c 7^ 0. With the new parameters 
and for A > we reduce to (14. 5p . 

The functions g of formula (14 .8p were found by searching for g such 
that the 4-dimensional system x = —g(x), y = —g'(x)y, where g G C 1 , 
g(0) = 0, (/(0) > 0, becomes superintegrable. More precisely, the 
system admits a first integral which is positive definite at the origin, 
so that it constrains the orbits near the origin on compact sets. That 
method only proves that all these functions g give a local isochronous 
center near the origin for 'x = —g(x). The problem of whether the 
center is indeed globally isochronous was not addressed in the original 
paper, but it can be also recovered by an analyticity argument, without 
appealing to Theorem 12.31 



5. DORIGNAC GLOBAL POTENTIAL REVISITED 

Let us consider 



(5.1) f(x,y) = (V 1 + 8e 3 ^' - l) (V 1 + 8e 3 ^ - l) -4, (x,y) 



G ^ , 



for the real parameters u>,(3 7^ 0. We have /(0,0) = 0, f(x,y) = 
f(y,x), and 

d 2 f(x, y) = 12/3c 3 ^ V Z—= - + 0, (x, y) G M 2 . 

VI + 8e 3 ^ 

The set f(x, y) = projects itself on the whole x-axis since the function 
x 1 — y \/T+ 8e 3l3x — 1 is a diffeomorphism of R onto (0, +00). By Theo- 
rem [37TJ, the equation f(x, y) = defines a global involution h : R — > R. 
Also now we can solve explicitly the equation and get 

1 2 eP x 
(5.2) h(x) = - In ; , x G R. 
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Next, formula (12.41) of Theorem 12.31 gives the following isochronous 
potential which is one of Dorignac's explicit global examples (see for- 
mula (35) in [7]) 
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